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随机 植物 短 纤 维 复 合 材料 珊 面 性 能 对 有 效 蛋 量 
和 拉 伸 行为 的 影响 * 


沈 琵 孙 晓 翔 刘 洋 
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摘 要 研究 了 界面 性 能 对 随机 短 云 杉 纤维 增强 聚 丙烯 (PP) 复 合 材料 宏观 拉 伸 性 能 的 影响 。 采 用 双 线 性 内 聚 力 模型 (CZMD) 描 
述 随 机 短 云 杉 纤维 和 PP 基体 间 非 理想 界面 的 力学 行为 , 建立 了 含 非 理 想 界面 的 随机 短 纤 维 增强 复合 材料 代表 性 单元 (RVE) 
的 二 维 有 限 元 模型 , 考虑 了 纤维 含量 、 长 细 比 、 随 机 分 布 和 随机 各 向 异性 弹性 以 及 PP 基体 弹 塑 性 的 影响 ; 模拟 了 不 同 纤维 含 
量 复 合 材料 的 实验 拉 伸 应 力 应 变 曲线 。 结 果 表 明 , 短 云 杉 纤维 /PP 基体 间 非 理想 界面 刚度 与 复合 材料 有 效 弹 性 模 量 之 间 有 
单调 递增 的 曲线 关系 , 即 E-K 曲线 ; 同一 复合 材料 不 同 纤 维 含 量 的 E-K 曲线 艇 有 一 个 临界 交点 。 在 交点 右 侧 强 界面 刚度 区 
复合 材料 有 效 模 量 随 着 纤维 含量 的 增加 而 提高 , 在 交点 左 侧 弱 界面 刚度 区 有 效 模 量 随 着 纤维 含量 的 增加 而 减少 。 三 种 不 同 
体积 含量 10%、20% 和 49% 的 云 杉 /PP 复合 材料 的 非 理 想 界面 刚度 可 用 E-K 曲线 和 实验 测 得 的 宏观 有 效 弹 性 模 量 确 定 , 云 杉 / 
PP 界面 初始 破坏 位 移 和 界面 完全 破坏 位 移 也 可 根据 模拟 拉 伸 应 力 应 变 曲线 确定 。 数 值 分 析 结 果 能 用 非 理 想 界面 刚度 来 解 
释 和 理解 随机 短 植物 纤维 体积 含量 对 复合 材料 宏观 有 效 模 量 的 影响 。 
关键 词 短 植物 纤维 , 拉 伸 性 能 , 有 限 元 , 纤维 /基体 界面 , 有 效 弹 性 模 量 , 内 聚 力 模型 
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ABSTRACT The influence of interface performance on the macroscopic tensile properties for the ran- 
dom short spruce fibers reinforced polypropylene (PP) composite materials was investigated. The me- 
chanical behavior of the imperfect interface between spruce fiber and PP matrix is described by the bilin- 
ear cohesive zone model (CZM), while a two-dimensional finite element model of the representative vol- 
ume element (RVE) with CZM for the material was developed in terms of the volume content, aspect ratio 
(AR) and random anisotropic elastic of random distribution short spruce fiber, as well as the influence of 
elastic plastic PP matrix. Experimental tensile stress strain curves for the composites with different fiber 
volume contents were simulated. The results show that there exists a common trend of monotone in- 
creasing for the curves of imperfect interfacial stiffness versus the effective modulus, namely E-K curves. 
The E-K curves for the composites with different volume fraction of fibers converged to a unique critical 
point (CP). In the range of higher interface stiffness the effective modulus of composites increase with the 
increase of fiber volume content, in the range of lower interface stiffness that is the opposite. For three 
spruce/PP composites with different fiber contents of 10%, 20% and 49% (in volume fraction), their imper- 
fect interfacial stiffness could be estimated by their E-K curves and the measured macroscopic effective 
elastic modulus through experiment. The displacement corresponding to the initial separation and that to 
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the entire separation of the interface could also be determined by the simulating tensile experimental 
stress strain curve of spruce/PP. Therefore, the results of numerical analysis base on the imperfect inter- 
facial stiffness can be used to explain and further understand the influence of random short fiber volume 
content on the effective modulus of spruce/PP composites. 

KEY WORDS plant short fiber, tensile behavior, finite element analysis, fiber/matrix interface, effective 


modulus, cohesive zone model 


植物 短 纤 维 增强 复合 材料 有 质 轻 、 环 保 和 可 
循环 使 用 等 优点 , 广泛 应 用 在 汽车 、 高 铁 列 车 和 
民用 飞机 等 领域 , 是 绿色 环保 复合 材料 的 研发 
热点 之 一 "”。 随 机 短 云 杉 纤维 增强 聚 丙 烯 复合 材 
料 ( 云 杉 /PP) 可 循环 使 用 , 是 具有 潜在 应 用 价值 的 
植物 纤维 增强 复合 材料 。 文 献 [3, 4] 研究 了 纤维 
微 结构 对 其 宏观 拉 伸 和 断裂 行为 的 影响 , 文献 [5] 
研究 了 该 材料 经 历 5 次 和 10 次 循环 使 用 后 的 拉 伸 
和 断裂 力学 性 能 。 植 物 短 纤 维 增强 复 合 材料 的 力 
学 性 能 主要 受 三 方面 因素 影响 , 即 纤维 和 基体 材 
料 的 力学 性 能 、 纤 维 的 微 结构 (尺寸 .长 细 比 、 体 
具 分 数 、 随 机 分 布 情况 等 ) 以 及 纤维 /基体 界面 的 性 


学 性 能 的 影响 。 

Hashin 提出 了 复合 材料 纤维 和 基体 间 的 非 理 想 
界面 的 概念 3, 考虑 该 界面 上 物理 量 如 应 力 和 位 移 的 
不 连续 , 来 研究 对 复合 材料 宏观 力学 性 能 影响 。 目 
前 , 描述 复合 材料 纤维 基体 非 理想 界面 力学 行为 的 模 
型 有 线性 弹 得 模型 (Linear spring model, LSM) 中 ,三 
相 材 料 模型 外 和 内 育 力 模型 (Cohesive zone model， 
CZM)" 09。CZM 模型 用 牵引 -分 离 定律 (Traction-Sep- 
aration law, T-S) 曲 线 描述 纤维 /基体 的 非 理想 界面 力 
学 行为 , 模拟 界面 硬化 、 软 化 直至 破坏 的 过 程 。 在 理 
论 和 数值 模型 中 工 S 曲线 被 广 为 应 用 。 目 前 还 没有 


能 。 其 中 纤维 /基体 界面 在 纤维 与 基体 之 间 传 递 
应 力 , 在 复合 材料 中 有 极为 重要 的 作用 "。 植 物 纤 
维 具 有 较 强 的 吸水 性 而 聚合 物 基体 具有 斥 水 性 ， 
影响 复合 材料 纤维 与 基体 之 间 的 厅 合 程度 出 现 
纤维 基体 界面 脱 粘 、 孔 洞 等 缺陷 , 使 载荷 难以 在 基 
体 与 纤维 之 间 进 行 有 效 传递 , 影响 复合 材料 的 宏 
观 力学 性 能 。 因 此 , 为 了 云 杉 /PP 复合 材料 的 优化 
设计 应 该 研究 其 界面 力学 行为 对 复合 材料 宏观 力 
学 性 能 的 影响 。 

一 些 学 者 研究 植物 短 纤 维 增强 复合 材料 力学 
性 能 时 发 现 , 材料 的 纤维 体积 含量 超过 某 一 数值 后 
其 力学 性 能 下 降 ”"。Hu 在 实验 中 发 现 "% 纤维 体积 
含量 为 50% 的 大 麻 纤 维 增强 聚 乳酸 复合 材料 的 强 
度 和 模 量 均 低 于 40% 材 料 的 相应 值 。 用 扫描 电子 
显微镜 观察 材料 的 断裂 表面 , 发 现在 纤维 体积 含量 
为 50% 的 材料 中 一 些 纤维 与 基体 并 没有 很 好 地 结 
合 。Arib 在 菠萝 叶 增 强 聚 丙烯 复合 材料 的 拉 伸 实 
验 中 发 现 ", 当 纤 维 体积 含量 小 于 10.8% 时 材料 的 
宏观 弹性 模 量 随 纤维 体积 含量 的 增加 而 上 升 , 而 当 
纤维 体积 含量 为 16.2% 时 其 宏观 弹性 模 量 反而 降 


有 效 和 通用 的 实验 方法 测定 CZM 界面 模型 参数 , 特 
别 是 对 随机 短 纤维 增强 复合 材料 的 界面 。 喻 明 等 

试 凑 法 估计 了 球形 颗粒 增强 复合 材料 界面 力学 
CZM 参数" 通过 有 限 元 模拟 研究 材料 的 破坏 和 界 
面 的 损伤 。Sreeranganathan 建立 三 维 颗粒 增强 复合 
材料 有 限 元 模型 , 研究 了 增强 相 和 基体 材料 的 力学 
性 能 以 及 材料 微 结 构 对 其 宏观 力学 性 能 的 影响 于 
但 是 没有 考虑 界面 因素 。 
使 用 内 聚 力 模型 描述 云 杉 /PP 复合 材料 的 非 理 
想 界 面 力 学 行为 , 已 通过 实验 与 数值 混合 法 初步 探 
索 了 该 随机 短 纤 维 复合 材料 宏 微 观 内 聚 力 参数 的 测 
量 方 法 "%; 用 两 个 细 观 理论 模型 探讨 了 随机 短 纤 维 
(颗粒 ) 增 强 复 合 材 料 的 非 理 想 界面 刚度 和 其 有 效 模 
量 的 理论 估计 方法 中 。 本 文采 用 数值 方法 , 基于 内 
聚 力 模型 表征 纤维 /基体 非 理 想 界 面 力 学 行为 , 进 

步 建立 具有 不 同 纤维 含量 的 随机 短 云 杉 /PP 复合 材 
料 的 代表 性 单元 (RVE) 的 有 限 元 模型 , 并 考虑 云 杉 
纤维 各 癌 异 性 弹性 性 能 、 纤 维 随 机 材料 主 方向 和 随 
机 位 置 , PP 基体 为 弹 塑 性 材料 。 研 究 纤 维基 体 间 界 
面相 性 能 对 随机 短 云 杉 /PP 复合 材料 宏观 拉 伸 行 关 
的 影响 。 使 用 对 比 实验 曲线 的 方法 得 到 相关 CZM 


低 。Hargitai 研究 大 麻 增 强 聚 丙烯 复合 材料 的 力学 
性 能 发 现 "", 当 材 料 的 体积 分 数 大 于 50% 时 材料 的 
弹性 模 量 出 现下 降 的 趋势 。 这 说 明 , 对 于 高 纤维 体 
积 含量 材料 界面 的 作用 极为 重要 。 植 物 纤维 复合 
材料 的 界面 性 能 的 相关 参数 , 特别 是 对 短 植 物 纤 
维 , 不 易 在 实验 中 直接 测量 和 表征 , 所 以 需要 建立 
考虑 纤维 基体 界面 因素 的 复合 材料 力学 模型 , 用 数 
值 方法 定量 分 析 以 研究 界面 性 能 对 其 材料 宏观 力 


Th 


性 


参数 , 并 定量 地 分 析 讨论 界面 刚度 与 材料 有 效 弹 
模 量 的 关系 。 


ce 
| 


1 模型 描述 
1.1 短 云 杉 /PP 复合 材料 
本 文 研究 的 几 种 随机 短 云 杉 纤维 增强 聚 丙烯 复 
合 材 料 ( 云 杉 /PP), 云 杉 纤 维 体积 含量 分 别 为 60%、 
49%、20% 和 10%。 云 杉 纤 维 是 用 工业 废料 (锯末 ) 加 


一 个 一 


0 


202303.10684v1 


chinaXiv 


aN 
| 


ChinaXiv 合 作 期 刊 


工 而 成 , 其 直径 和 长 度 的 统计 平均 尺寸 分 别 约 为 
0.4 mm 和 0.8 mm, 长 细 比 4R 约 为 2。 聚 丙烯 (PP) 基 
体 材料 的 型 号 为 APPRYL3400 MA1。 根 据 云 杉 纤 
维 的 形态 特点 外 可 将 其 视 为 正 交 各 向 异性 弹性 材 
料 , 加 工 方式 使 得 其 材料 主 方向 一 般 与 纤维 轴线 方 
向 的 夹 角 随 机 分 布 ; PP 基体 是 均匀 各 向 同性 弹 塑 性 
材料 。 用 奥 林 帕 斯 工业 显微镜 得 到 的 云 杉 /PP 复合 
材料 的 表面 细 观 形 貌 如 图 1(a) 所 示 。 本 文 用 到 的 云 
杉 和 PP 材料 的 弹性 常数 列 于 表 125 0 PP 的 实验 应 
力 应 变 曲 线 在 文献 [3] 中 给 出 。 用 Instron 1343 电 液 伺 
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1 云 杉 /PP 复合 材料 的 表面 细 观 形 貌 和 RVE 模型 的 
几何 及 边界 条 件 

Fig.1 Microstructure surface image (a) and geometry and 
boundary conditions of RVE model (b) of spruce/ 
PP composites 


表 1 云 杉 纤维 和 PP 基体 的 弹性 模 量 


Table 1 Elastic properties of spruce fiber and PP matrix 


Properties Symbol ee PP 
fiber 
E 10991.0 
E/ MPa E, 716.0 1573.0 
E; 435.0 
Vu 0.42 
V Vi 0.48 0.4 
Va 0.50 
Gy 724.0 
G/ MPa G's 557.0 
Gs 31.6 
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服 静 动态 试验 机 进行 云 杉 /PP 复合 材料 的 拉 伸 试验 , 使 
标准 板 型 拉 伸 试 件 , 试验 段 截面 尺寸 为 10x4 mm， 
在 室温 下 用 位 移 控制 加 载 (速率 为 20 mm/min), 测 出 
云 杉 纤维 的 体积 含量 分 别 为 49%、20% 和 10% 的 复 
合 材 料 的 应 力 应 变 曲线 。 
1.2 考虑 界面 的 RVE 有 限 元 模型 

针对 云 杉 /PP 复合 材料 的 特点 , 本 文 建立 了 可 以 
考虑 云 杉 纤维 体积 含量 各 向 异性 弹性 性 能 、 纤 维 随 
机 材料 主 方向 和 随机 位 置 等 微 结 构 信 息 , 上 且 PP 基 体 
为 弹 塑性 材料 的 随机 短 云 杉 /PP 复合 材料 代表 性 单 
元 (RVE) 的 有 限 元 数值 模型 。 

模型 中 假设 : 1) 随 机 短 云 杉 纤 维 是 均匀 各 向 异 
性 弹性 材料 , 在 基体 材料 内 的 分 布 位 置 、 分 布 方向 是 
随机 的 , 其 材料 主 方向 一 般 与 纤维 轴线 方向 是 随机 
的 ; 2) 基 体 是 均匀 各 向 同性 弹 塑性 材料 ; 3) 云 杉 /PP 
间 界 面 是 没有 脱 粘 的 非 理 想 界面 , 其 力学 行为 用 
CZM 界面 模型 来 描述 ; 4)RVE 内 含有 足够 多 的 短 纤 
维 , 纤维 、 基 体 的 体积 含量 满足 复合 材料 的 混合 定 
律 ; 5) 破 坏 均 发 生 在 纤维 /基体 界面 处 , 而 基体 和 纤 
维 均 不 发 生 破坏 。 为 简化 起 见 , 设 受 力 前 该 RVE 中 
各 处 、 所 有 纤维 的 纤维 /基体 界面 性 能 皆 相 同 , 即 所 
有 界面 的 CZM 模型 参数 相同 。 

随机 短 云 杉 /PP 复合 材料 的 二 维 RVE 模 型 如 图 
1 所 示 , 为 使 模型 的 纤维 微 结构 特征 接近 真实 的 复 
合 材 料 , 纤维 分 布 位 置 和 纤维 轴线 取向 是 随机 的 , 用 
数字 1.2 表示 每 个 云 杉 纤 维 在 xy 平面 内 的 两 个 材料 
主 方向 , 如 图 1 所 示 , 建 模 时 材料 主 方向 1 与 纤维 轴 
线 的 夹 角 也 是 随机 的 。 每 个 云 杉 纤 维 通过 周围 的 一 
层 以 CZM 表征 的 接近 零 厚度 非 理想 界面 与 PP 基体 
粘 结 。 
图 1 中 模型 的 长 与 宽 分 别 用 a 和 b 表 示 , 在 本 文 
中 RVE 模 型 尺寸 参照 Lam 等 外 的 云 杉 /PP 复合 材料 
拉 伸 试验 试 件 的 均匀 实验 段 尺 寸 , 即 二 维 RVE 模型 
厚度 为 4 mm, 取 a=8 mm, b=5 mm, 纤维 长 细 比 为 
2。 在 单 向 拉 伸 时 , RVE 模 型 的 边界 条 件 可 以 表示 为 


=- bb 
u(0,y)=0 ye| .| 


u(a,y)= Au yel- 如 旨 (1 
0,0=0 
其 中 Az 为 拉 伸 量 , 本 文 取 Aw=0.3 mm。 

为 生成 RVE 有 限 元 模型 以 及 界面 相 单元 的 前 
处 理 , 本 文 用 随机 序列 吸附 (Random Sequential Ab- 
sorption, RSA) 方 法 辣 , 用 Fortran 语言 编写 前 处 理 程 
序 , Cohesive 界面 单元 采用 与 纤维 和 基体 单元 共 结 
点 的 方式 植 入 RVE 有 限 元 模型 。RSA 方 法 的 基本 
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思想 是 : 随机 取向 的 纤维 随机 分 布 于 基体 之 中 , 新 生 
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在 本 文 的 RVE 有 限 元 模型 中 , 采用 ABAQUS 程 


成 的 纤维 不 得 与 其 它 已 生成 的 纤维 产生 谍 入 或 接触 
的 情况 , 否则 去 除 已 生成 纤维 并 重新 生成 新 的 纤维 ， 
直至 材料 当前 的 纤维 体积 分 数 达 到 目标 体积 分 数 。 
在 本 文中 , 前 处 理 程序 实现 RSA 的 方法 是 , 每 当 新 
生成 纤维 单元 时 其 占用 的 有 限 元 网 格 节点 不 得 被 已 
生成 的 纤维 占用 , 否则 重新 赋予 纤维 的 位 置 和 取向 ， 
继续 判断 该 纤维 是 否 符合 要 求 , 直至 生成 符合 要 求 
的 纤维 , 代 入 在 相应 的 位 置 上 。 
1.3 纤维 基体 界面 的 Cohesive 单 元 表征 
本 文 将 随机 短 云 杉 /PP 复合 材料 的 纤维 /基体 界 
面 看 作 非 理想 界面 , 用 内 聚 力 模型 (CZM) 表 征 , 非 理 
想 界 面 的 内 聚 力 t 和 界面 分 离 量 6 间 的 关系 用 CZM 
的 牵引 -分 离 定律 , 即 Ts 曲线 来 表征 。 
目前 细 观 力学 实验 方法 如 纤维 拔 出 实验 、 微 滴 
脱 粘 实验 . 单 丝 段 裂 实验 和 压 头 顶 出 实验 等 能 测量 
界面 的 某 些 力学 性 能 , 但 是 对 于 纤维 长 细 比 小 的 随 
机 短 云 杉 /PP 复合 材料 , 还 没有 适合 的 实验 方法 测量 
真实 的 纤维 /基体 界面 内 聚 力 模型 的 工 S 曲线 。 为 了 
突出 界面 力学 行为 的 主要 特征 、 减 少 分 析 参 数 , 本 文 
CZM 模型 使 用 双 线 性 的 工 S 曲线 , 如 图 2 所 示 , 其 中 
界面 分 离 量 6 小 于 初始 破坏 位 移 6 时 为 弹性 阶段 
(OA 段 ), 表示 Cohesive 界面 刚度 , 最 大 牵引 力 如 为 
Cohesive 界面 强度 ; 8 为 界面 完全 破坏 位 移 , 当 6,< 
6< 6 时 , 界面 进入 软化 阶段 (AB 段 ), 界面 发 生 损 伤 ， 
牵引 力 1 逐渐 下 降 , 6 二 6 时 , 界面 达到 完全 脱 粘 的 状 
态 外 (BC 段 )。 在 TS 曲线 的 弹性 阶段 , 二 维 t 与 6 的 
本 构 关系 为 

i 

| 


I 
式 中 下 角 标 na 和 s 分 别 表 示 界 面 的 法 线 和 切线 方向 ,各 
和 分 别 表 示 界 面 法 向 和 切 向 的 牵引 力 , 2 和 & 分 别 
表示 界面 法 向 和 切身 的 分 离 位 移 , Ks 和民, 分别 表 示 
法 向 和 切 向 的 界面 刚度 , 为 Cohesive 界面 厚度 。 
t/MPa K: Interfacial stiffness 


下: Softening Interfacial stiffness 
tnax: Maximum traction 


G: Fracture energy 
6o: Initial separation 
6r: Debonding separation 
OA: Elastic (stage1) 
AB: Solftening (stage2) 
BC: Debonding (stage3) 


0 60 6f 6/mm 
Separation 


2 cohesive 单 元 双 线 性 牵引 -分 离 曲 线 
Fig.2 Bilinear traction-separation law of cohesive element 


序 中 的 二 维 4 节 点 Cohesive 单元 (COH2D4) 模 拟 纤 
维 / 基 体 非 理想 界面 的 力学 行为 。 对 于 双 线 性 CZM 
模型 的 TS 曲线 , 可 用 三 个 独立 参数 为 一 组 确定 其 
牵引 -分 离 关 系 , 即 界面 刚度 KK、 相 应 的 初始 破坏 位 
移 6 和 和 完全 破坏 位 移 5.。 

为 简化 起 见 , 计算 中 假设 Cohesive 界面 的 法 向 
和 切 向 的 TS 曲线 相同 , 相应 各 个 双 线 性 CZM 参数 
相同 , 即 K=K,=K.、6.=6%w=6, 和 6=6i=6 等 等 。 

试 算 和 验证 表明 , Cohesive 界面 (单元 ) 厚 度 对 
复合 材料 拉 伸 应 力 应 变 曲 线 的 模拟 影响 极 小 故 可 以 
忽略 , 计算 中 马 值 取 一 个 很 小 的 值 0.5 mm。 由 于 界 
面相 所 占 的 体积 很 小 , 在 计算 纤维 含量 时 不 考虑 界 
面相 体积 的 影响 。 

本 文采 用 上 述 RVE 有 限 元 模型 , 模拟 云 杉 /PP 
复合 材料 的 拉 伸 实验 应 力 -应 变 曲线 , 反复 迭代 试 
算 , 确定 其 纤维 /基体 界面 内 聚 力 模型 的 一 组 参数 ， 
即 及 .和 6 等 。 

1.4 RVE 有 限 元 模型 计算 的 实施 及 验证 

使 用 ABAQUS6.10 有 限 元 分 析 软 件 计 算 分 析 
随机 短 云 杉 /PP 复合 材料 RVE 二 维 有 限 元 数值 模 
型 , 模拟 云 杉 纤维 和 PP 基体 的 有 限 单元 错 采 用 4 节 
点 平面 应 力 缩减 积分 单元 (CPS4R), 界面 单元 使 用 4 
节点 Cohesive 单元 (COH2D4)。 由 于 RVE 有 限 元 模 
型 内 含有 足够 多 的 随机 分 布 且 取 向 随机 的 短 纤维 ， 
有 限 元 网 格 划 分 必须 足够 密 , 才能 得 到 精确 的 数值 
结果 。 为 了 保证 弹 塑 性 、 非 线性 计算 的 收敛 , 且 有 足 
够 的 精度 , 必须 验证 RVE 有 限 元 网 格 密度 对 计算 得 
到 的 复合 材料 宏观 有 效 性 能 , 如 对 宏观 有 效 模 量 预 
测 的 数值 结果 影响 。 

对 纤维 含量 为 V=20% 的 RVE 有 限 元 模型 , 图 3 
给 出 了 四 种 CZM 界面 刚度 KK, 分 别 为 5x10; MPa/mm、 


1620 
上 40x24 、 80*48 160x96 320x192 480x384 
上 一 v v v 
1570F 40x24 -80%48 160x96 320x192 480x384 
上 全 一 全 一 全 人 A 全 
1520 上 上 
区 —— K=5x10? MPa/mm 
&_ 1470 上 本 
三 上 —®— A=7x10 MPa/mm 
ww 1420F 一 全 _ K=4x10* MPa/mm 
上 —¥— K=10° MPa/mm 
1370 上 
F 40*24 80x48 160x96 320x192 480x384 
1320 5 ee A lb 
1270 FF 40*24 80*48 160x96 ”320x192 480x384 
[上 是 加 | 
1220 C5 Ps 1 i a | 1 下 汪汪 
10” 103 104 105 105 


Number of elements 


图 3 网 格 密度 对 预测 的 宏观 有 效 模 量 影响 
Fig.3 Effective modulus of models with different mesh 
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7x10 MPa/mm、 4x10+ MPa/mm、 和 1 x10° MPa/mm, 
同时 考虑 五 种 不 同 的 有 限 元 网 格 密度 , 即 分 别 是 
40x24、80x48、160x96、320x192 和 480x384, 计算 得 到 
的 宏观 有 效 模 量 。 试 算 结 果 表明 , 当 网 格 密度 大 于 80x 
48 时 宏观 有 效 模 量 预测 的 数值 结果 偏差 小 于 0.1%。 
此 , 本 文选 择 适中 的 单元 网 格 密度 为 160x96, 纤维 和 
基体 单元 对 应 的 尺寸 为 0.05 mmx0.052 mm, 以 保证 
计算 精度 。 

上 述 试 算 结 果 还 表明 , 由 于 云 杉 /PP 复合 材料 的 
纤维 长 细 比 很 小 (AR= 2), 且 纤 维 是 各 向 异性 的 , 材 
料 主 方向 是 随机 的 , 短 纤 维 轴 向 的 取向 随机 性 影响 
很 小 。 在 不 影响 预测 结果 的 情况 下 , 为 了 减少 有 限 
元 网 格 划 分 的 困难 和 避免 单元 网 格 过 密 , 建 模 时 短 
纤维 轴线 的 取向 只 随机 分 布 在 试 件 的 轴 向 (x 方向 ) 
和 其 垂直 方向 , 如 图 4 所 示 。 图 4 给 出 的 是 纤维 体积 
含量 分 别 为 10%、20% 和 49% 的 随机 短 云 杉 纤维 增 
强 聚 丙烯 复合 材料 RVE 的 有 限 元 模型 网 格 和 CZM 
界面 单元 划分 图 。 在 纤维 体积 含量 49% 的 RVE 模 
型 中 , 基体 和 纤维 单元 总 数 为 15360 个 , cohesive 单 
元 数 为 2832 个 , 总 节点 数 为 18685 个 。 

按照 上 述 含有 界面 单元 的 RVE 有 限 元 方法 建 
立 的 随机 短 云 杉 /PP 复合 材料 的 “数值 试 件 ”即使 纤 
维 ( 及 其 含量 )、 基 体 和 界面 的 力学 性 能 相同 , 由 于 每 
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其 中 和 分别 为 有 效 弹 性 模 量 的 统计 平均 值 和 方 
差 , 即 j 为 材料 在 该 界面 刚度 下 的 数值 计算 得 到 的 
有 效 弹 性 模 量 E。 


2 结果 和 讨论 


本 文 基于 双 线 性 内 聚 力 模型 表征 的 云 杉 /PP 复 
合 材 料 的 纤维 基体 界面 性 能 , 包括 非 理想 界面 刚度 、 
界面 强度 和 界面 断裂 能 , 与 之 对 应 的 的 双 线 性 CZM 
模型 的 独立 参数 分 别 是 非 理想 界面 刚度 KKK、 初始 破 
坏 位 移 & 和 完全 破坏 位 移 8。 
2.1 非 理 想 界面 刚度 与 有 效 模 量 的 单调 递增 关系 

为 了 研究 随机 短 云 杉 /PP 复合 材料 中 纤维 基体 
间 的 非 理想 界面 刚度 KK 对 其 宏观 有 效 弹 性 模 量 E 的 
影响 , 先 利用 本 文 建立 的 随机 短 云 杉 /PP 复合 材料 考 
虑 纤维 基体 界面 行为 的 代表 性 单元 (RVE) 有 限 元 数 
值 模 型 , 通过 对 界面 刚度 天 的 参数 研究 ,用 FEM 数 
值 模型 得 到 宏观 有 效 弹 性 模 量 5 与 非 理想 界面 刚度 
KK 的 关系 曲线 。 

分 别 用 10%、20% 和 49% 三 种 纤维 体积 含量 的 
云 杉 /PP 材料 的 RVE 有 限 元 数值 模型 , 模拟 相应 体 
只 含量 真实 材料 的 拉 伸 力学 行为 , 其 中 云 杉 纤维 的 


个 “数值 试 件 ” 中 短 纤维 的 位 置 、 取 向 和 材料 各 向 异 
性 主 方向 分 布 是 随机 的 , 不 同 的 数值 试 件 预测 得 到 
的 复合 材料 宏观 有 效 模 量 的 数值 是 不 同 的 ,有 一 定 
的 分 散 性 。 参 照 文献 [23] 的 方法 , 如 在 预测 同一 纤 
维 含量 的 云 杉 /PP 复合 材料 宏观 有 效 模 量 时 , 给 定 同 
一 个 界面 刚度 和 其 它 材料 参数 后 , 本 文 用 10 个 短 纤 
维 的 位 置 . 取 向 和 材料 各 向 异性 主 方向 随机 分 布 的 
不 同 “ 数 值 试 件 ” 预 测 宏观 有 效 模 量 , 并 用 高 斯 概率 


密度 函数 中 进行 统计 分 析 , 即 
(a) ©® 


Interface 
Matrix 
4 代表 体 RVE 模 型 的 网 格 划分 


Fig.4 Mesh of RVE model with AR=2: (a) V=10%, (b) 
VYF20% and (c) V=49% 


长 细 比 AR 为 2, 纤维 材料 为 随机 正 交 各 向 异性 弹 
性 ; 选取 3x10; MPa/mm~10” MPamm 之 间 共 20 个 
不 同 界 面 刚度 值 , 计算 每 个 界面 刚度 值 对 应 的 10 
个 复合 材料 “数值 试 件 ” 的 有 效 弹 性 模 量 , 得 到 对 应 
上 述 每 个 界面 刚度 值 的 统计 平均 有 效 弹 性 模 量 羽 
并 以 PP 基体 的 弹性 模 量 ,为 参照 , 得 到 无 量 纲 化 
的 有 效 模 量 BE, 和 非 理 想 界面 刚度 KK 之 间 的 关系 , 即 
云 杉 /PP 复合 材料 的 已 开关 系 曲线 (E- 玉 曲线 )。 上 述 
三 种 纤维 体积 含量 的 云 杉 /PP 材料 BE-K 曲线 , 在 图 5 
中 给 出 。 


~ Critical Point, CP 
Experimental value 


WwW 0.8 C IF10% 
ee oO Vz20% 
0.6 性 49% 
Imperfect interface Perfectinterface 
/=10% -oO- VF10% 
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K/MPa-mm'! 
5 云 杉 /PP 复合 材料 的 BE-K 曲线 及 各 含量 的 实验 值 


Fig.S E-K curves of spruce/PP composite and experimen- 
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图 5 表明 , 对 于 同一 体积 含量 的 云 杉 /PP, 其 有 效 
模 量 E 会 随 着 其 非 理 想 界面 刚度 K 值 的 增加 而 单调 
增加 。 对 于 不 同体 积 含量 的 该 复合 材料 , 可 得 到 一 
簇 E-K 曲 线 。 云 杉 /PP 的 E-K 关 系 曲线 , 可 分 为 两 个 
阶段 : (D) 当 界面 刚度 玉 值 较 大 , 如 图 中 当 界 面 刚度 
增加 到 10 MPa/mm, 或 趋向 无 穷 时 , 复合 材料 界面 
接近 理想 粘 接 情况 , 各 含量 的 B-K 曲 线 均 趋 于 各 自 
的 水 平 渐进 线 ; (2) 界 面 刚度 KK 值 较 小 时 , 曲线 中 部 
上 升 阶段 是 典型 的 非 理 想 界 面 , 描述 了 云 杉 /PP 复合 
材料 在 非 理想 界面 刚度 影响 下 的 宏观 有 效 模 量变 化 
趋势 , 即 非 理想 界面 刚度 与 宏观 有 效 模 量 之 间 具 有 
单调 递增 的 关系 。 

文献 [3] 实 验 测 量 了 10%、20% 和 49% 三 种 纤维 
体积 含量 的 云 杉 /PP 材料 的 应 力 应 变 曲 线 , 可 以 得 到 
其 宏观 有 效 模 量 EE,。 在 此 利用 本 文 含 界 面 的 RVE 
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E 随 着 体积 含量 Vi 的 增加 而 增加 是 有 条 件 的 , 与 其 
纤维 /基体 的 非 理想 界面 刚度 值 的 大 小 有 关 : 

(1) 在 CP 的 左下 方 , 非 理 想 界面 刚度 值 小 于 临 
界 界 面 刚度 K,, 有 效 模 量 巨 低 于 基体 弹性 模 量 瓦 ， 
即 提高 短 纤维 含量 未 起 到 增强 的 效果 , 有 效 模 量 了 
随 着 纤维 体积 模 量 Vi 的 增加 而 减 小 , 此 区 域 为 弱 界 
面 刚度 区 ; 

(2) 各 体积 含量 的 曲线 汇 交 于 CP, CP 的 纵 坐 标 
值 等 于 基体 弹性 模 量 5,, 即 无 量 纲 化 的 纵 坐 标 值 等 
于 证 

(3) 在 CP 右上 方 , 非 理想 界面 刚度 值 大 于 临界 
界面 刚度 及, 有 效 模 量 E 高 于 基体 弹性 模 量 E,, 有 
效 模 量 E 随 着 短 纤维 体积 含量 Vi 的 增加 而 增加 , 此 
区 域 为 强 界面 刚度 区 。 对 于 云 杉 /PP 材料 , 当 界 面 刚 
度 增 加 到 10 MPa/mm, 各 含量 的 B-K 曲线 均 趋 于 各 


且 


+ 


有 限 元 模型 得 到 的 云 杉 /PP 复合 材料 的 E-K 曲线 , 进 


面 刚度 KK, 如 表 2 所 示 。 实 验 测 得 的 有 效 模 量 和 相 
应 界面 刚度 对 应 的 E-K 曲 线 上 点 的 坐标 为 (K, Ew/ 
万 ,), 在 图 5 中 分 别 用 三 个 空心 点 来 表示 。 

在 本 文 的 基本 假设 条 件 下 , 因 非 理想 界面 刚度 
与 宏观 有 效 模 量 之 间 有 单调 递增 的 关系 , 用 这 个 数 
值 方 法 通过 实验 测 得 的 随机 短 纤维 复合 材料 的 宏观 
有 效 模 量 Es 以 估计 其 纤维 /基体 间 非 理想 界面 刚度 
尺 的 值 。 
2.2 纤维 含量 和 界面 刚度 对 有 效 弹 性 模 量 的 影响 

从 图 5 可 以 看 出 , 同一 短 纤维 复合 材料 不 同 纤 
维 含 量 的 E-K 曲 线 徐 名 汇 交 于 一 点 。 称 这 个 交点 为 
纤维 含量 影响 复合 材料 有 效 模 量 的 临界 交点 (Criti- 
cal Point CP), 把 临界 交点 对 应 的 复合 材料 有 效 模 
量 称 为 临界 有 效 模 量 E,, 把 该 点 对 应 的 纤维 /基体 非 
理想 界面 刚度 称 为 临界 界面 刚度 Ks,。 从 图 5 可 知 , 对 
于 纤维 长 细 比 为 2 的 云 杉 /PP 复合 材料 , 其 临界 有 效 模 
量 玉 =1594 MPa, 临界 界面 刚度 Ks=107 MPa/mm。 

临界 交点 CP 表明 , 随机 短 纤维 增强 复合 材料 的 
纤维 /基体 非 理 想 界 面 刚度 影响 纤维 体积 含量 Vi 对 
宏观 有 效 模 量 E 的 增强 或 减弱 趋势 。 即 5-K 关 系 曲 
线 簇 交点 CP 表明 , 随机 短 纤维 复合 材料 的 有 效 模 量 

表 2 Cohesive 参 数 及 材料 的 有 效 弹性 模 量 实验 值 

Table 2 Cohesive parameters and experimental effective 
elastic modulus 


VYr/%  E:/MPa K/(MPa/mm) /mm Oo/mm 


10 1341.4 2600 0.0220 0.93 
20 1201.3 3960 0.0105 1:10 


49 1432.2 25360 0.0016 0.80 


自 的 水 平 渐进 线 , 此 时 相当 于 纤维 与 基体 界面 是 理 
想 粘 接 的 情况 。 
对 于 本 文 三 种 纤维 体积 含量 为 10% 、20% 和 
49% 的 云 杉 /PP 材料 , 由 表 2 中 实验 测 得 的 宏观 有 效 
模 量 ,估计 的 纤维 /基体 间 非 理想 界面 刚度 KK 的 
值 , 均 小 于 临界 界面 刚度 K,, 位 于 临界 交点 左 侧 的 
弱 界 面 刚度 区 。 纤 维 体积 含量 10% 和 20% 云 杉 /PP 
材料 的 非 理想 界面 刚度 KK 值 非常 接近 , 用 上 述 对 已 
KK 曲线 簇 交 点 的 分 析 , 如 图 5 所 示 , 可 以 解释 纤维 体 
积 含量 较 高 的 20% 云 杉 /PP 的 有 效 模 量 反 而 会 低 于 
10% 的 有 效 模 量 , 导致 材料 有 效 模 量 降低 的 原因 是 
弱 界 面 刚度 。 
因此 , 在 这 类 随机 短 纤 维 云 杉 /PP 复合 材料 的 以 
提高 有 效 模 量 为 目标 的 设计 中 , 除了 提高 纤维 体积 
含量 以 外 , 还 需要 特别 重视 纤维 /基体 的 界面 设计 ， 
特别 是 要 保证 非 理想 界面 刚度 值 位 于 临界 交点 CP 
右 侧 的 强 界面 刚度 区 。 
针对 各 向 同性 颗粒 增强 复合 材料 , 文献 [20] 采 
细 观 力学 的 Mori-Tanaka(M-T) 方 法 和 稀释 解 方法 
研究 了 非 理想 界面 刚度 对 该 复合 材料 有 效 模 量 的 影 
响 , 得 到 不 同 颗 粒 体积 含量 的 单调 递增 的 EK 关系 
曲线 艇 , 并 用 这 两 个 理论 模型 证 明了 E-K 关 系 曲 线 
秘 存 在 唯一 的 交点 。 本 文 的 数值 模型 扩展 了 上 述 
细 观 力学 模型 的 能 力 , 考虑 到 短 纤维 是 各 向 异性 材 
料 、 且 纤维 长 细 比 AR 大 于 1, 得 到 随机 短 纤维 增强 
复合 材料 5-K 关 系 曲 线 秘 和 临界 交点 , 以 便 基于 实 
验 测 得 的 宏观 有 效 模 量 预测 纤维 /基体 的 非 理 想 界 
面 刚度 。 
2.3 界面 CZM 参数 对 应 力 应 变 曲 线 的 影响 及 确定 
用 本 文 建立 的 随机 短 云 杉 纤 /PP 复合 材料 的 考 
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数值 模型 模拟 的 三 种 不 同 纤维 体积 含量 的 云 杉 /PP 


足够 多 的 随机 分 布 短 纤维 , 在 纤维 体积 含量 、 纤 维 尺 


寸 、 纤 维 长 细 比 AR 


况 下 , 模拟 该 材料 的 实验 拉 伸 应 力 应 变 曲线 通过 


对 数值 模拟 的 分 析 ， 


的 双 线 性 CZM 界面 模型 的 三 个 独立 参数 , 即 界面 网 


皆 与 复合 材料 实际 “一 致 ”的 情 


可 以 研究 纤维 基体 非 理 想 界面 


度 玉 ,界面 初始 破坏 位 移 & 和 界面 完全 破坏 位 移 


对 该 复合 材料 拉 伸 行为 的 影响 。 
在 应 力 应 变 曲线 的 线 弹 性 阶段 , 计算 结果 表明 ， 


在 弹性 阶段 CZM 模 型 的 6 和 6 参数 , 对 复合 材料 的 
宏观 弹性 模 量 巨 没有 影响 ; 复合 材料 的 纤维 基体 间 


CZM 的 非 理想 界 理 


刚度 天 只 影响 弹性 阶段 应 力 应 


变 曲线 的 斜率 , 即 复合 材料 的 宏观 弹性 模 量 已 上 
述 数值 分 析 表 明 , 随机 纤维 /基体 间 的 非 理想 界面 刚 
度 KK 与 该 复合 材料 的 宏观 弹性 模 量 E 之 间 有 单调 递 
增 关 系 。 因 此 , 由 实验 得 到 的 茶 一 纤维 体积 含量 的 


宏观 有 效 弹 性 模 量 ， 


可 根据 本 文 数值 模型 计算 得 到 


的 该 复合 材料 的 EB-K 关 系 曲线 唯一 地 确定 该 体积 含 


量 的 复合 材料 的 纤 


维 /基体 的 非 理 想 界 于 


刚度 玉 。 


用 该 模型 模拟 拉 伸 应 力 应 变 曲 线 , 弹性 阶段 可 与 


文献 [3] 中 得 到 的 实 


验 曲线 吻合 。 图 6 给 出 了 数值 


模拟 得 到 的 纤维 体积 含量 为 10% 的 云 杉 /PP 应 力 应 
变 曲线 , 用 上 述 方 法 再 得 到 了 非 理想 界面 刚度 到 


如 表 2 所 示 。 图 6 


FP 尽管 几 


数值 模 曲 线 的 CZM 


界面 的 6 和 6 值 不 同 , 但 其 弹性 阶段 多 与 实验 曲线 吻 


合 很 好 。 


在 应 力 应 变 曲线 的 非 线 弹性 阶段 , 确定 天 之 后 ， 


数值 结果 表明 , 如 民 


CZM 模型 参数 m 和 5 只 对 复合 材料 拉 伸 应 力 应 变 曲 
线 的 非 线 弹性 阶段 有 影响 , 即 曲线 屈服 和 强化 阶段 
的 强度 和 破坏 , 非 理 想 界 再 


6 所 示 纤 维 /基体 的 非 理 想 界 面 


刚度 玉 对 此 阶段 无 影响 。 


通过 模拟 应 力 应 变 曲线 的 非 线 弹性 阶段 , 可 以 


明 , 在 及 和 8 不 变 时 ， 


通过 数值 研究 初步 
含量 为 10% 的 云 杉 /PP 的 实验 应 力 应 变 曲线 , 结果 表 


角 定 参数 6 和 6:。 模 拟 纤 维 体积 


随 6 的 增加 , 拉 伸 曲线 最 大 应 力 


增加 (图 6a), 以 通过 对 实验 曲线 的 模拟 确定 该 材料 


的 CZM 界面 初始 破坏 位 移 6。 
天 和 上 不 变 时 , 随 着 8 的 变化 拉 人 


试 算 结果 还 表明 , 在 
H 线 强化 段 的 形状 


随 之 变化 (图 6b), 可 通过 对 实验 曲线 的 模拟 确定 该 


材料 的 CZM 界面 完全 破坏 位 移 &。 从 而 得 到 纤维 


体积 含量 为 10% 的 云 杉 /PP 复合 材料 的 界面 CZM 参 


数 6 和 6 的 值 , 如 表 2 所 示 。 


同 理 , 根据 云 杉 /PP 的 纤维 体积 含量 为 20% 和 


49% 的 实验 应 力 应 变 全 曲线 中 


j 本 文 方法 确定 相应 


的 纤维 /基体 界面 的 CZM 参数 , 如 表 2 所 示 。 用 本 文 


复合 材料 的 拉 伸 曲线 与 实验 曲线 的 比较 , 如 图 7 所 
示 。 由 有 限 元 模拟 与 实验 曲 


线 比 较 而 确定 的 三 种 体 


站 含量 的 云 杉 /PP 复合 材料 
面 刚度 玉 、 初 始 破坏 位 移 & 


| 界面 的 CZM 参数 , 即 界 
和 完全 破坏 位 移 6 的 值 , 


如 表 2 所 示 。 使 用 本 文 建立 的 随机 短 云 杉 /PP 复合 
材料 的 考虑 纤维 基体 非 理想 界面 的 RVE 有 限 元 模 


型 , 可 模拟 该 材料 的 不 同 的 纤维 体积 含量 实验 拉 人 


应 力 应 变 曲 线 。 可 以 看 出 ， 


型 后 有 限 元 的 预测 精度 提高 , 数值 预测 与 文献 [3] 的 


引入 非 理 想 界面 CZM 模 


实验 结果 吻合 。 
45 
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6 Cohesive 参数 对 材料 应 力 应 变 


线 的 影响 


Fig.6 Influence of cohesive parameters to stress- Strain 
curves (Vr=10%) (a) 6, (b) 6 
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图 7 应 力 应 变 曲线 数值 与 实验 结果 


Fig.7 Stress- strain curves achieved by numerical and ex- 


perimental method 
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微 结构 (尺寸 、 长 细 比 、 体 积 含量 
特点 , 建立 了 该 
限 元 模型 ,采用 双 线 性 内 聚 力 模型 描 
材料 的 非 理 想 界 二 
向 异性 弹性 性 


3 结 论 


1. 本 文 基于 随机 短 云 杉 /PP 复合 材料 的 纤维 的 
` 随 机 分 布 情况 等 ) 
复合 材料 的 代表 性 单元 (RVE) 的 有 
述 云 杉 /PP 复合 
行为 , 同时 考虑 了 云 杉 纤维 的 各 
能 、 纤 维 随机 材料 主 方向 和 随机 位 置 ， 


HH 


PP 基体 为 弹 塑性 材料 。 


2. 用 含有 非 理想 界面 的 RVE 有 限 元 模型 通过 


对 云 杉 /PP 的 界面 刚度 的 参数 研究 结果 表明 , 随机 短 
云 杉 /PP 复合 材料 的 有 效 模 量 巨 与 云 杉 /PP 间 非 理想 


界面 刚度 天 值 有 单调 递增 的 曲线 关系 , 即 EF-K 曲 
线 。 不 同 纤 维 体积 含量 的 云 杉 /PP 复合 材料 的 E-K 


| 线 艇 有 一 个 临界 交点 (CP), CP 从 纤维 /基体 的 非 理 


4D 


ai 


俱 


想 界 面 刚度 大 小 的 角度 , 反映 了 纤维 体积 含 
合 材 料 宏观 有 效 模 量 影响 的 两 种 趋势 , 即 在 CP 点 左 
1 弱 界 面 刚度 区 , KK 值 偏 小 ,有效 模 量 EE 随 着 纤维 体 


只 模 量 WW 的 增加 而 减 小 ; 在 CP 点 右 侧 强 界面 刚度 区 


量 对 复 


KK 值 偏 大 , 有 效 模 量 E 随 着 纤维 体积 模 量 VW 的 增加 
而 增加 , 当 界 面 刚度 天 值 无 限 大 时 界面 为 理想 粘 接 ， 


短 
的 变化 可 以 用 现 有 的 到 


可 用 实 


纤维 复合 材料 的 有 效 模 量 E 随 着 纤维 体积 模 量 多 
E 论 公式 计算 "。 
3. 对 于 某 一 纤维 体积 含量 的 云 杉 /PP 复合 材料 ， 
验 得 到 的 有 效 拉 伸 弹性 模 量 , 由 本 文 数值 模 


去 


型 得 到 的 E-K 曲线 , 预测 该 短 纤维 /基体 间 非 理想 界 


面 的 刚度 KK 值 。 


云 杉 /PP 复合 材料 , 在 模拟 其 实验 拉 人 


4. 针对 三 种 不 同体 积 含量 10%、20% 和 49% 的 


应 力 应 变 曲 线 


基础 上 , 使 用 本 文 的 RVE 数值 模型 和 方法 可 分 别 得 
到 它们 的 EK 曲线 , 初步 确定 了 云 杉 /PP 基体 间 
CZM 模 型 的 参数 , 即 非 理 想 界 面 刚 度 K, 界面 初始 


破坏 位 移 6, 和 界 凋 
想 界面 刚度 和 E-K 1 


完全 破坏 位 移 6.。 用 预测 的 非 理 
| 线 , 解释 了 纤维 体积 


含量 较 高 


的 20% 云 杉 /PP 的 有 效 模 量 反 而 会 低 于 10% 的 有 效 
模 量 的 实验 现象 。 
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